A Hamiltonian Krylov-Schur-type method based on the symplectic Lanczos process
نویسندگان
چکیده
We discuss a Krylov-Schur like restarting technique applied within the symplectic Lanczos algorithm for the Hamiltonian eigenvalue problem. This allows to easily implement a purging and locking strategy in order to improve the convergence properties of the symplectic Lanczos algorithm. The Krylov-Schur-like restarting is based on the SR algorithm. Some ingredients of the latter need to be adapted to the structure of the symplectic Lanczos recursion. We demonstrate the efficiency of the new method for several Hamiltonian eigenproblems.
منابع مشابه
Solving Large-scale Quadratic Eigenvalue Problems with Hamiltonian Eigenstructure Using a Structure-preserving Krylov Subspace Method
We consider the numerical solution of quadratic eigenproblems with spectra that exhibit Hamiltonian symmetry. We propose to solve such problems by applying a Krylov-Schur-type method based on the symplectic Lanczos process to a structured linearization of the quadratic matrix polynomial. In order to compute interior eigenvalues, we discuss several shift-and-invert operators with Hamiltonian str...
متن کاملA Krylov-schur Approach to the Truncated Svd
Computing a small number of singular values is required in many practical applications and it is therefore desirable to have efficient and robust methods that can generate such truncated singular value decompositions. A new method based on the Lanczos bidiagonalization and the Krylov-Schur method is presented. It is shown how deflation strategies can be easily implemented in this method and pos...
متن کاملKrylov subspace methods for the Dirac equation
The Lanczos algorithm is evaluated for solving the time-independent as well as the time-dependent Dirac equation with arbitrary electromagnetic fields. We demonstrate that the Lanczos algorithm can yield very precise eigenenergies and allows very precise time propagation of relativistic wave packets. The unboundedness of the Dirac Hamiltonian does not hinder the applicability of the Lanczos alg...
متن کاملPerturbation Analysis of Hamiltonian Schur and Block-Schur Forms
In this paper we present a complete perturbation analysis for the Hamiltonian Schurform of a Hamiltonian matrix under similarity transformations with unitary symplectic matrices.Both linear asymptotic and non-linear perturbation bounds are presented. The same analysis isalso carried out for two less condensed block-Schur forms. It suggests that the block forms areless sensitive ...
متن کاملLanczos-type Solvers for Non-hermitian Linear Systems
In this overview we discuss iterative methods for solving large linear systems with sparse (or, possibly, structured) nonsymmetric (or, non-Hermitian) matrix that are based on the Lanczos process. They feature short recurrences for the generation of the Krylov space and for the sequence of approximations to the solution. This means low cost and low memory requirement. For very large sparse non-...
متن کامل